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Abstract. We present optical high-resolution spectra of a 
sample of sixteen post-AGB stars and IRC -I-I0216. Of the 
post-AGB stars, ten show C2 Phillips (Ainu-Xil]+) and 
Swan (d^Bg - a^B^ ) and CN Red System (A^B - X'^T.+) 
absorption, one CH"*" (A^B — X^S"*") emission, one CB+ 
absorption, and four without any molecules. We find typ- 
ically Trot 43 - 399, 155 - 202, and 18 - 50 K, logiV - 
14.90 - 15.57, 14.35, and 15.03 - 16.47 cm'^ for C2, CB+, 
and CN respectively, and 0.6 < N (CN) / N {C2) < 11.2. 
We did not detect isotopic lines, which places a lower 
limit on the isotope ratio of ^^C/^^C > 20. The pres- 
ence of C2 and CN absorption is correlated with cold dust 
(Jdust < 300 K) and the presence of CH"*" with hot dust 
(Idust > 300 K). All objects with the unidentified 21 /im 
emission feature exhibit C2 and CN absorption, but not 
all objects with C2 and CN detections exhibit a 21 /im 
feature. The derived expansion velocity, ranging from 5 to 
44 km s"""^, is the same as that derived from CO millimeter 
line emission. This unambiguously proves that these lines 
are of circumstellar origin and are formed in the AGB 
ejecta (circumstellar shell expelled during the preceding 
AGB phase). Furthermore there seems to be a relation 
between the C2 molecular column density and the expan- 
sion velocity, which is attributed to the fact that a higher 
carbon abundance of the dust leads to a more efficient 
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acceleration of the AGB wind. Using simple assumptions 
for the location of the molecular lines and molecular abun- 
dances, mass-loss rates have been derived from the molec- 
ular absorption lines and are comparable to those obtained 
from CO emission lines and the infrared excess. 

Key words: molecular processes — circumstellar matter 
- stars: AGB and post-AGB — line: identification 



1. Introduction 

The first study of molecules in the optical spectrum of a 
post-AGB star concerned the presence of C3 absorption 
and C2 emission in the reflected light of the lobes of the 
Cygnus Egg Nebula (Crampton et al. 1975). Renewed in- 
terest was triggered by the discovery by Waelkens et al. 
(1992), Balm & Jura (1992), and Hall et al. (1992) of 
CH+ emission in the optical spectrum of the famous Red 
Rectangle. Recently, Waelkens et al. (1995) reported on 
the presence of CH+ absorption in HD 213985. In this pa- 
per (which we will refer to as Paper II) we will study the 
presence of C2, CN, and CH+ in the optical spectrum of 
thirteen (post-)AGB stars. 

Bakker et al. (1996c, Paper I) were the first to analyze 
the molecular bands of C2 and CN in the spectrum of 
the post-AGB star HD 56126 and showed that the expan- 
sion velocity and the excitation conditions are consistent 
with the lines being formed in the AGB ejecta. Extend- 
ing this work, Bakker et al. (1995) showed that the same 
is valid for the molecular absorption lines in four other 
post-AGB stars (IRAS 04296-1-3429, IRAS 05113-^1347, 
IRAS 08005-2356, and AFGL 2688) and argued on the 
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basis of the relation between CO (or OH for IRAS 08005- 

2356) and C2 or CN expansion velocity that these molecu- 
lar absorption lines are formed in the AGB ejecta and are 
therefore of circumstellar origin. Hrivnak (1995) discussed 
the presence of C2 and C3 in the low-resolution spectra of 
nine post-AGB stars. 

Theoretical models of the formation and dissociation 
of molecules in the extended envelope of the carbon-rich 
AGB star IRC -M0216 by Cherchneff et al. (1993) showed 
that, e.g., C2 and CN are only present in a thin shell 
of material within the extended envelope. Close to the 
star carbon and nitrogen are locked up in complex stable 
molecules such as C2H2 and HCN, while at larger dis- 
tances the interstellar ultraviolet radiation field photodis- 
sociates these molecules to C2 and CN. At even larger 
distances the UV radiation field photodissociates simple 
molecules into their constituent atoms and ions. The net 
effect is that simple molecules exist only in a thin shell of 
material. 

Interstellar C2 has been discussed extensively by van 
Dishoeck & Black (1982). They have shown that the ex- 
citation of C2 is a sensitive balance between photoexcita- 
tion and coUisional (de-)excitation. By modeling the exci- 
tation, the relative population over the rotational energy 
levels of the C2 ground state can be used to determine the 
particle density, radiation field, and kinetic temperature 
of the line-forming region. In a separate paper (Paper III 
in preparation) we will apply this model to the stars stud- 
ied, while an earlier account of this work can be found in 
Bakker et al. (1995). 

In § 2 we discuss the criteria used in selecting our sam- 
ple of mainly carbon-rich post-AGB stars, and describe 
the observations and data reduction. § 3 describes the 
method used in identifying molecular bands and determi- 
nation of the expansion velocities of the AGB ejecta. Ro- 
tational diagrams are used to derive molecular rotational 
temperatures, column densities, and mass-loss rates. The 
results are discussed (§ 4) for the sample as a whole and 
for each star separately. We will finish with a short con- 
clusion (§5). 

2. Observations and data reduction 

2.1. Selection of program stars 

After the serendipitous discovery of C2 and CN in the op- 
tical spectrum of HD 56126 (Paper I) we started an exten- 
sive observing campaign in order to find additional post- 
AGB stars with molecular absorption or emission lines. 
Assuming that the presence of carbon-based molecules 
is related to the high carbon abundance of the AGB 
ejecta, we have selected those stars with a carbon-rich 
circumstellar environment. The presence of the 3.3 and 
3.4-3.5 /xm Polycyclic Aromatic Hydrocarbon (PAH) fea- 
tures and of the unidentified 21 /zm feature (Kwok et al. 
1989) were used as criteria for the carbon-rich nature of 



the AGB ejecta (Table 1). Recently the number of ob- 
jects exhibiting the 21 /im feature has been extended 
(e.g., Henning et al. 1996, Justtanont et al. 1996), and 
these objects have been selected for a follow-up study. The 
sample was supplemented with those stars of special im- 
portance for the theory of post-AGB evolution (e.g., the 
metal-depleted post-AGB binaries: HD 52961, HR 4049, 
BD +39°4926, Red Rectangle, and HD 213985). Two O- 
rich post-AGB stars (e.g, HD 161796 and IRAS 08005- 
2356) were added to see whether 0-rich star have the 
carbon based molecules, and the well studied carbon 
star IRC-l-10216 was added to the list, since its spec- 
trum could possibly be used as a template for iden- 
tifying molecular features. With a limiting magnitude 
of rriv = 14 this resulted in a list of sixteen Post- 
AGB stars (and IRC -1-10216) observable from La Palma. 
IRAS 05341-1-0852 was added to the list at a later stage. 
The optical spectrum and the molecular bands are de- 
scribed by Reddy et al. (1997). 

The observed stars are listed in Table 1, and the re- 
sulting spectra in three selected windows are presented 
in Figs. 1,2, and 3. To facilitate comparison all velocities 
given are heliocentric (marked with the symbol Vq) and 
the corrections were obtained with the Starlink RV utility 
Version 2.2. 

2.2. Observations and data reduction 
2.2.1. WHT/UES 

The observations were made in four different runs between 

February 1992 and August 1994 with the Utrecht Echelle 
Spectrograph (UES) mounted on the Nasmyth platform 
of the 4.2m WiUiam Herschel Telescope (WHT) on La 
Palma. The echellograms of February 1992 were recorded 
on an EEV CCD-05-30 detector with 1242 x 1152 pixels of 
22.5 X 22.5 /um^ each, while in 1994 a Tektronix TK1024A 
device was used with 1024 x 1024 pixels of 24 x 24 ^m^ 
each. The wavelength coverage of the UES ranges from 
3000 to 11000 A, and even using the echelle with 31.6 fines 
per mm, one needs several settings because of the small 
frame size of the CCD's. For our observations we used 
settings with central wavelengths of 4020, 5261 and 7127 
A, giving a wavelength coverage from about 3650 to 10000 
A, with some overlap. For each setting Tungsten flatfield 
and Thorium-Argon (Th-Ar) calibration exposures were 
taken. 

The echellograms from the UES were reduced by TS 
using the echelle package of IRAF V2.9 running on a 
DEC3100 workstation. All frames were bias corrected us- 
ing the overscan columns of the CCD's. The traced orders 
were extracted with the optimal extraction method using a 
variance weighting algorithm. Subsequently, the extracted 
spectra were divided by the identically extracted fiatfields 
and normalized to the continuum. The wavelength cali- 
bration was performed using the FIGARO list of Th-Ar 
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Table 1. Infrared fluxes and spectral features of the seventeen program stars. 



Object 


Sp.T. 


my 


/12 [Jy] 


/25 [Jy] 


3.3 * 


3.4-3.5 * 


21 * 


Remark 


IRAS 04296+3429 


GO la 


14.21 


12.74 


45.94 


e 


e 


e 




IRAS 05113+1347 


G8 la 


14.40 


3.78 


15.30 


e 


nd 


c 




IRAS 05341+0852 


F4 I 


12.8 


4.51 


9.85 


e 


e 


c 


s-process enhanced 


HD 44179 


B9 I 


8.84 


421.60 


456.10 


e 


nd 


nd 


IRAS 06176-1036; Red Rectangle 


HD 52961 


F8 I 


8.50 


4.53 


2.22 


no 


no 


nd 


IRAS 07008+1050 


HD 56126 


F5 I 


8.23 


24.51 


116.70 


c 


nd 


c 


IRAS 07134+1005 


IRAS 08005-2356 


F5 1 


11.46 


17.96 


51.80 


no 


no 


nd 


OH maser 


IRC +10216 


C9.5 


11.00 


47530.00 


23070.00 


no 


no 


ud 


IRAS 09452+1330 (carbon star) 


HR 4049 


Fl I 


5.52 


48.25 


9.55 


e 


nd 


nd 


IRAS 10158-2844 


HD 161796 


F3 lb 


7.27 


6.12 


183.50 


no 


no 


nd 


IRAS 17436+5003 (oxygen rich) 


IRAS 20000+3239 


G8 la 


13.40 


15.03 


70.97 


no 


no 


c 




AFGL 2688 


F5 lae 


14.00 


339.00 


3041.00 


e 


no 


e 


Egg Nebulae 


IRAS 22223+4327 


GO la 


13.30 


2.12 


37.10 


no 


e 


e 




HD 235858 


G2 la 


9.30 


73.88 


302.40 


6 


e 


6 


IRAS 22272+5435 


HD 213985 


B9 I 


8.83 


5.57 


4.66 


te 


nd 


nd 


IRAS 22327-1731 


ED +39°4926 


Fl I 


9.24 


nd 


nd 


no 


no 


no 




IRAS 23304+6147 


G2 la 


13.15 


11.36 


59.07 


no 


no 


e 





*: e: emission; nd: not detected; no: not observed, 3.3, 3.4-3.5, and 21 fj,m infrared features 



lines. About 700 calibration lines, evenly distributed over 
the orders, were used in a two-dimensional fit with typical 
rms residuals of 0.25 km s^^ (about 0.1 pixel). By checking 
the wavelength of a number of telluric oxygen absorption 
lines, the drift between calibration and target frames was 
found to be less than 1 km s~^. 

The resolution was determined from a number of iden- 
tified telluric lines (Moore et al. 1966) yielding a FWHM 
of the line profile of 6.0 ± 1.0 km and a spectral re- 
solving power of i? ^ 5 X 10^. The typical signal-to-noise 
ratio is SNR = 100, and the minimum equivalent width 
detectable is about 7 mA (see App. A only at CDS). 

2.2.2. McDonald/CS21 

In search for the CH"*" band towards HD 56126 (no such 
band was detected), we have obtained high-resolution 
spectra using the 2DC0UDE (CS21) spectrograph (Tull 
et al. 1995) of the 2.7m telescope at the McDonald obser- 
vatory. Light was fed to the TK3 CCD, 2048 x 2048 pixels 
with each 24 x 24 /um^, at the Fl focus using E2 echelle 
having 52.6759 grooves mm"-'^. The spectra were reduced 
by EJB using the data reduction package IRAF: bias and 
scattered light subtracted, flat fielded and wavelength cal- 
ibrated using a ThAr arc spectrum. The wavelength cal- 
ibration has a rms internal error of 0.42 mA which cor- 
responds to 0.016 km at the wavelength of interest. 
The resolution was determined from the Thorium lines as 
160,000. 

3. Molecular absorption and emission lines 



Table 2. Circumstellar molecular bands in the optical spectra 

of p<)st-:VGB stars. 



Molecule 


System 


Band 




{v',v") 


A ran 

[A] 


ge 


12^12^ 


Phillips 


A^n, 




(1,0) 


10137 


10213 


12^13^ 


Phillips 


A^n, 




(1,0) 


10166 


10242 


12^12^ 


Phillips 


A^Hu 


-X^S+ 


(2,0) 


8753 


8913 


12(213(2 


Phillips 


A^Hu 


-XiE+ 


(2,0) 


8797 


8957 


12(-,12(-, 


Phillips 


A^Hu 




(3,0) 


7717 


7809 


12(213(2 


Phillips 


A^n, 


-xis+ 


(3,0) 


7769 


7861 


12(212^ 


Phillips 


A^Hu 




(4,0) 


7910 


7002 


12(213(2 


Phillips 


A^Hu 


-X^EJ 


(4,0) 


7965 


7057 


12(212(2 


Swan 


d^H, 


-a^Hu 


(0,0) 


5100 


5166 


12^13^ 


Swan 


d^Hg 


-a^Hu 


(0,0) 


5995 


5161 


12^12^ 


Swan 




-a^Hu 


(1,0) 


4710 


4740 


12(213(2 


Swan 


d^Hg 


-a^Hu 


(1,0) 


4705 


4735 


12(2l4j^ 


Red System 


A^H- 


-x2e+ 


(1,0) 


9130 


9220 


«C"N 


Red System 


A^n- 


-x2e+ 


(1,0) 


9159 


9249 


12^14^ 


Red System 


A^H- 




(2,0) 


7874 


7918 


13^14^ 


Red System 


A^H- 


-x2e+ 


(2,0) 


7918 


7962 


12^14^^ 


Red System 


A^H- 




(3,0) 


6927 


6954 


13^14^ 


Red System 


A^H- 


-x2e+ 


(3,0) 


6967 


6994 


12^14^^ 


Red System 


A^H- 




(4,0) 


6190 


6225 


13(214^ 


Red System 


A^H- 


-x2e+ 


(4,0) 


6244 


6279 


12(2ljj+ 




A^H- 


-XiE+ 


(0,0) 


4220 


4280 


i3(2iH-h 




A^H- 


-XiE+ 


(0,0) 


4220 


4280 


i2(2iH+ 




A^H- 


-XiE+ 


(1,0) 


3950 


3990 


13clH+ 




A^H- 


-XiE+ 


(1,0) 


3950 


3990 
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7720 



7740 



7760 
A [A] 



7780 



7800 



Fig. 1. Spectra of the observed stars in the wavelength region of the C2 A^IIu — X^S^ (3,0) band. The top spectrum is a 
synthetic spectrum computed using Trot = 200 K, log AT = 15.60 cm~^ and ii ~ 5 x 10^. All other spectra are corrected for 
the velocity of the molecular lines such that the molecular lines are at their rest wavelengths. The three absorption lines near 
7775 A arc due to the OI(l) multiplet. Note that the wavelength part from 7740 A to 7770 A is not observed for most of the 
objects in the sample. 



3.1. Identification of molecular hands 

Molecular absorption line profiles resemble closely those 
of telluric absorption lines, while molecular emission line 
profiles can easily be confused with cosmic spikes. To pre- 
vent confusion with cosmic spikes most integrations were 
divided into two or more shorter exposures, and the indi- 
vidual spectra were combined after the reduction process 
using a median filter to remove cosmic spikes. A number 
of (hot) stars have been observed to identify telluric lines. 

In this paper we have restricted ourselves to twelve 
molecular bands and their -"^"^C isotopic species (Table 2) 
of three different species (C2, CN, and CH+). Table 3 
gives for each star observed and for each molecular 
band, whether it was detected (in absorption or emis- 
sion), tentatively detected (mostly when the lines are 
severely blended by photospheric lines), or not detected 
at all. We have carefully studied all spectra to find ad- 



ditional "prominent" molecular bands, but none were 
detected. Of course this does not exclude the presence 
of weak molecular bands with equivalent widths of less 
than 10 ml. IRAS 05113-M347, IRAS 20000-^3239, and 
HD 235858 show also CN Red System band photospheric 
absorption (Fig. 2). These lines are easily distinguisable 
from the circumstellar CN lines, since the lines are much 
broader and they are identified with much higher energy 
levels (typically 30 < A^" < 90). 

We have carefully looked for the isotopic lines, e.g., 
i^C^^C, and 13CH+, with negative results. From 

this we deduce a typical lower limit on the isotope ratio 
of 12C/13C > 20. 

In Paper I we pointed out that optical depth effects 
play an important role for the stronger molecular bands; 
therefore in this study we have limited the analysis to the 
weaker molecular bands: C2 PhiUips (3,0) (Fig. 1), CN 
Red System (3,0) (Fig. 2), and CH+ (0,0) (Fig. 3). 
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Fig. 2. Spectra of the observed stars in the wavelength region of the CN A^II — X'^S^ (3,0) band. The top spectrum is a 
synthetic spectrum of circumstellar ^^CN computed using Trot = 25 K, logN = 16.50 cm~^ and J? ~ 5 x 10"', followed by a 
synthetic spectrum of photospheric ^^CN using Trot = 5500 K and log = 16.78 cm~^. All spectra are corrected for the velocity 
of the molecular lines such that the circumstellar molecular lines arc at their rest wavelengths. Most of the unidentified narrow 
absorption lines (see HR 4049 which has no circumstellar or photospheric lines in this part of the spectrum) are due to telluric 
O2 and H2O (Moore et al. 1966) or photospheric CN (only for HD 235858, IRAS 20000+3239, and IRAS 05113+13477). 
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Fig. 3. Spectra of the observed stars in the wavelength region of the CH"'" A^II — X^E"*" (0,0) band. The top spectrum is a 
synthetic spectrum computed using Trot = 200 K, log AT = 14.70 cm~^ and R = 5 x 10^. All other spectra are corrected for 
the velocity of the molecular lines such that the molecular lines are at their rest wavelengths. Note the broad photospheric 
absorption features in the spectrum of HD 213985 (e.g., Srll(l) and the somewhat narrower FeII(27) feature.) 



The wavenumbers of the C2 Phillips bands have been 
taken from Chauville et al. (1977) ((1,0) and (2,0)) and 
Ballik & Ramsay (1963) (3,0) and CH+ (0,0) and (1,0) 
bands from Carrington & Ramsay (1982). Conversion 
from wavenumber in vacuum to wavelength in air was 
made by by applying Cauchy's formula. Wavelengths of 
the CN Red System were extracted from the SCAN tape 
(J0rgenson & Larsson 1990). 



3.2. Determination of the expansion velocity, rotational 
temperature, column density, time-scale, and mass- 
loss rate 

3.2.1. Expansion velocity 

The first step in analyzing a molecular band is to check 
that all the candidate molecular lines yield the same radial 
velocity. The average radial velocity of all lines is a good 
measure for the true velocity of the line-forming region. 



For absorption lines this is the line-of-sight to the star, 
while for emission lines it can be anywhere within the 
"slit" (defined as the area on the CCD used for extracting 
the spectrum) of the telescope. 

In Paper I we showed that the molecular absorption 
lines in the spectrum of HD 56126 are formed in the AGB 
cjccta. Now we investigate for the current sample if there 
is a correlation between the expansion velocities of the 
AGB ejecta derived from CO (or OH for IRAS 08005- 
2356) emission line profiles, and the velocity differences 
between the molecular absorption lines and the system 
velocities. The system velocity was taken to be the cen- 
tral velocity of the CO or OH line emission profile or from 
radial velocity studies on binaries (Table 4). The veloci- 
ties of the photospheric lines are not a good estimate of 
the system velocity, since most of the star in the sample 
are somewhat pulsating and/or are binaries. 

The upper panel of Fig. 4 shows that the expansion ve- 
locity derived from CO or OH emission lines and those de- 
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Table 3. The detection of molecular bands in the optical spectra of post-AGB stars. 



Object 


C2 S 


wan 


C2 


Phillips* 




CN Red System 


* 


CH 


+ * 






"^-i-i- 
-a^Hu 


A^n 






A^n- 






i 1 -I— r 

A^H - 






(0,0) 


(1,0) 


(1,0) 


(2,0) 


(3,0) 


(1,0) 


(2,0) 


(3,0) 


(4,0) 


(0,0) 


(1,0) 


IRAS 04296+3429 


no 


no 


a 


a 


a 


a 


a 


a 


b 


no 


no 


IRAS 05113+1347 


a 


b 


a 


a 


a 


a 


a 


a 


b 


no 


no 


IRAS 05341+0852 


no 


no 


a 


a 


a 


a 


a 


a 


no 


no 


no 


HD 44179 


np 


np 


np 


np 


np 


np 


np 


np 


np 


e 


e 


HD 52961 


np 


np 


np 


np 


np 


np 


np 


np 


np 


np 


np 


HD 56126 


a 


a 


a 


a 


a 


a 


a 


a 


a 


np 


np 


IRAS 08005-2356 


ta 


no 


a 


a 


a 


a 


a 


np 


np 


b 


b 


IRC +10216 


b 


b 


a 


a 


a 


a 


ta 


b 


b 


no 


no 


HR 4049 


np 


np 


np 


np 


np 


np 


np 


np 


np 


np 


np 


HD 161796 


no 


no 


np 


np 


np 


np 


np 


np 


np 


np 


np 


IRAS 20000+3239 


a 


b 


no 


a 


a 


a 


a 


a 


a 


no 


no 


Ar KjLi zDOo 


a 


I- 

D 


e 


a 


a 


a 


a 


a 


a 


no 


no 


IRAS 22223+4327 


a 


a 


no 


a 


a 


a 


a 


a 


a 


no 


no 


HD 235858 


a 


a 


a 


a 


a 


a 


a 


b 


b 


no 


no 


HD 213985 


np 


np 


np 


np 


np 


np 


np 


np 


np 


a 


a 


BD +39°4926 


no 


no 


no 


no 


no 


no 


no 


no 


no 


np 


no 


IRAS 23304+6147 


a 


a 


no 


a 


a 


a 


a 


a 


a 


no 


no 



* a: absorption; b: blended; e: emission; no: not observed; np: not present; ta: tentative absorption 
C is '^C, N is "N, and H is 'H 



Table 4. Adopted system and AGB outflow velocities of programs stars. 



Object ''hn^.Im- ' 'Vy^. • ' 'Vx|. ' llt'iiiark lli'forouco 



IRAS 04296+3429 


-66.0 


-59.0 


12.0 


CO(J = 2-1) 





IRAS 05113+1347 


-12.0 


2.0 




photospheric lines 


B 


IRAS 05341+0852 


9.6 


25.0 




photospheric lines 


R 


HD 44179 


2.7 


20.7 


6.0 


binary P = 310 ± 3 days 


V,J 


HD 52961 








pulsator P = 70.8 days 


F 


HD 56126 


71.0 


85.6 


10.0 


CO 


Z 


IRAS 08005-2356 


46.7 


61.1 


49.8 


OH maser 


T 


IRC +10216 


-26.2 


-19.1 


14.1 


pulsator P = 638 days 


H,D 


HR 4049 


-44.5 


-32.9 




binary P = 429 days 


V 


HD 161796 


-36.0 


-0.4 


11.5 


CO(J = 1-0) 


L 


IRAS 20000+3239 


14.0 


-4.1 


12 


CO(J = 1-0) 





AFGL 2688 


-33.3 


-49.2 


22.8 




Y 


IRAS 22223+4327 


-30.0 


-42.2 


14.0 


C0(1 -0) 





HD 235858 


-30.9 


-43.1 


11.6 


CO 


Z 


HD 213985 


-42.5 


-45.7 




binary P = 259 days 


V 


BD +39°4926 








binary P = 775 days, no infrared excess 


K 


IRAS 23304+6147 


-15.9 


-25.8 


15.5 


C0(2 - 1) 


W 



Visr is the system velocity in the local standard of rest system [km s~^] 

Vq is the system velocity in the heliocentric system [km s~^] 
" Wexp the expansion velocity of the AGB ejecta [km s~^] 
B: this study; D: Dyck et al 1991; F: Fernie 1995; H: Huggins et al. 1988; K: Kodaira et al. 1970; J: Jura et al. 1995; L: Likkel 
ct al. 1987; O: Omont ot al. 1993; R: Rcddy ct al. 1997; T: Te Lintel Hekkert et al. 1991; V: van Winckel et al. 1995; W: 
Woodsworth et al. 1990; Y: Young et al. 1992; Z: Zuckerman et al. 1986 



Table 5. Physical parameters derived from the molecular lines. 









C2 A^Hu 


-X^E+ (3; 


0) or (2,0) ^ 












CN A^H - 


X^E+ (2,0) or 


(3,0) '> 










JV(CN) 
,^V(C2) 






Object 


Vq 


^exp 




logiV^oi 


logM 






Vq 


''cxp 


T,ot log 


A^mol " 


log M 




5v 


e 


12 






[km s 


ii 
J 


[km s"'] 


[K] 




[cm-^] 


[Mq yr-i] 


[km s 




[km s"'] 


[K] 


[cm-^] 


[Mq yr-i] 


















±2.0 






+0.10 


±1 










±2.0 




±0.10 


±1 








±2.5 






IRAS 04296+3429 


-66.7 ± 0. 


,4 


7.7 


138 + 14 


15 


.21(15.20) 


-5.8 


—63, 


.6 


+ 0.6 


4.6 


25 ±06 


15.91 


-5.2 


-3.1 


±0, 


.8 


4.9 


> 


20 


IRAS 05113+1347 


-2.5 + 0, 


,6 


4.5 


198 + 31 


15 


.13(15.13) 


-6.1 


-3, 


,6 


+ ] 


L.l 


5.6 


19 + 02 


15.89 


-5.1 


1.1 


+ 1, 


,3 


5.6 


> 


20 


IRAS 05341+0852 


13.4 + 0, 


,8 


11.6 


77 + 09 


15 


.00(14.61) 


-6.0 


15, 


.0 


+ 0.4 


10.0 


38 + 04 


15.86 


-5.0 


-1.6 


+ 1, 


.0 


7.2 






HD 56126 


77.3 + 0, 


,1 


8.3 


242 + 34 


15 


.31(15.28) 


-5.5 


78, 


.3 


+ 0.2 


7.3 


25 + 08 


15.51 


-5.2 


-1.0 


+ 0, 


.2 


1.5 


> 


20 


IRAS 08005-2356 


17.4 + 0, 


.4 


43.7 


149 + 10 


15 


.55(15.60) 


-4.7 


22, 


.3 


±0.2 


42.2 


38 + 02 


15.48 


-4.7 


-4.9 


±0, 


.4 


0.8 


> 


11 


IRC +10216 


-33.3 ± 0, 


,7 


14.2 


43 + 16 


14 


.90(15.15) 


-6.3 


-33, 


.8 


+ 0.5 


14.7 


28 ±05 


15.03 


-6.0 


-0.5 


±0, 


.9 


1.5 






IRAS 20000+3239 


-16.9 + 0, 


,4 


12.8 


234 + 33 


15 


.41(15.44) 


-5.3 


-17, 


,0 


+ 0.3 


12.9 


50 + 02 


16.47 


-4.1 


0.1 


+ 0, 


,5 


11.2 


> 


20 


AFGL 2688 


-66.5 ± 0, 


,3 


17.3 


56 + 08 


15 


.47(15.51) 


-4.8 


-64, 


.5 


+ 0.5 


15.3 


18 + 03 


16.10 


-4.1 


-2.0 


+ 0, 


.6 


4.2 


> 


19 


IRAS 22223+4327 


-57.2 ± 0, 


,2 


15.0 


399 + 36 


15 


.57(15.55) 


-4.6 


-55, 


.0 


±0.3 


12.8 


30 + 04 


16.34 


-3.8 


-2.2 


±0, 


.4 


5.8 


> 


20 


HD 235858 


-52.2 + 0, 


,3 


9.1 


119 + 35 


15 


.20(15.27) 


-5.7 


-51, 


.8 


±0.1 


8.7 


24 ±17 


15.03 


-5.8 


-0.4 


±0, 


.3 


0.6 


> 


11 


IRAS 23304+6147 


-39.7 + 0, 


,3 


13.9 


281 + 21 


15 


.57(15.57) 


-5.0 


-39, 


,2 


±0.5 


13.4 


43 ±14 


16.17 


-4.3 


-0.5 


±0, 


.6 


3.8 


> 


20 








CH+ 


A'n-X'E+(0,0) 
































HD 44179' 


17.3 + 0, 


,5 


3.4 


202 + 16 


14 


.44(14.49) 




























> 


22 


HD 213985 


-52.4 + 0, 


,6 


6.7 


155 + 22 


14, 


.35(14.35) 




























> 


38 



C2 (2,0) for IRC +10216, else C2 (3,0) 
^ CN (2,0) for IRAS 08005-2356, HD 235858, and IRC +10216, else CN (3,0) 

Vexp = Vsya — l'C2/CN/CH + 

in brackets: sum over J" or N" levels of observed transitions 

" Sv = VCN — VC2 

' the emission line spectrum has been analyzed with A^H J' = as energy zerolevel. 
Vexp and Trot are real, but the column density of HD 44179 is rather meaningless 
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Fig. 4. Top panels: expansion velocities derived from CO (or OH for IRAS 08005-2356) versus expansion velocities derived 
from C2 and CN. Lower panels: C2 and CN column density versus C2 and CN expansion velocity. The solid lines in the upper 
two panels give the relations in which the two velocities are equal. The dashed lines are linear least-squares fits as given in the 
text. IRAS 05341-1-0852 is not plotted. 



rived from C2 or CN absorption lines are identical within 
the estimated error of 2 km s~^. The expansion velocities 
determined from C2, CN, CH+ arc tabulated in Table 5. 
This proves unambiguously that the line-forming region 
of these molecular absorption lines is the AGB ejecta (cir- 
cumstellar shell). The optical molecular lines are not re- 
solved {FWHM ^ 6.0 km s~^), which puts an upper limit 
of 6.0 km s~^ on the turbulent broadening and velocity 
stratification within the line-forming region. The CH"*" de- 
tections are discussed in §4, but to first order they have the 
same range of expansion velocities. A linear least-squares 
fit (excluding IRAS 05113+1347, IRAS 05341-^0852 for 
which no CO data are available) gives: 



(CO) = 1.12 X Uexp (C2) + 0.36 (1) 
Ve^p (CO) = 1.10 X t;exp (CN) + 1.80 (2) 

with correlation coefficients of 0.98 and 0.97, respectively. 
All velocities are in units of km s~^. 



3.2.2. Rotational temperature and molecular column den- 
sity 

The second step is to determine the equivalent widths 
of the lines and to construct rotational diagrams (Fig. 5, 
6, 7, and 8). All rest wavelengths, (absorption) oscilla- 
tor strengths and equivalent widths of the lines used are 
tabulated in App. A (available only at CDS or from the 
authors). The adopted method is extensively discussed in 
Paper I, and we refer the interested reader to that paper 
for all details on the rotational diagram. Unfortunately 
there is a error in Eq. 5 of Paper I. The correct equation 
is: S% = (2 J" + 1) (page 208, Hcrzberg 1950). 

The (absorption) oscillator strengths were computed 
in the usual way (see Paper I for further references) us- 
ing band oscillator strengths of /(2.0) = 1-44 x 10~^ and 
/(3,o) = 6.672 X 10-'' for C2 Phillips (2,0) and (3,0) band 
respectively (Langhoff et al. 1990, Langhoff 1996), and 
/(o_o) = 5.45 X 10-3 for CH+ (0,0) (Larsson & Siegbahn 
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Fig. 5. Rotational diagrams for the stars showing the Phillips C2 (3,0) absorption. Note that for those stars for which transitions 
from high J" levels are observed (e.g., IRAS 22223+4327), the rotational diagram shows a flattening due to optical pumping of 
the molecule (see §3.2.2 for details). 



1983). For CH+ emission the emission oscillator strength 
is given by gj' f emission = 9 J"f absorption- AU data of the 
CN Red System were extracted from the SCAN tape of 
J0rgenson & Larsson (1990) with the oscillator strength 
multiplied by 0.734 as suggested by the authors. 

The rotational temperature is determined using rota- 
tional constant, By, values of 1.8111, 1.8907, 13.9302, and 
11.4227 cm-i for C2 v" = (Marenin & John- 

son 1970), CN X2E+ v" = (Brocklchurst ct al. 1971), 
CH+ XiS+ v" = (Carrington & Ramsay 1982) and 
CH+ A^Iiv' = (Carrington & Ramsay 1982), respec- 
tively. Higher order corrections terms (e.g. D„, H^) where 
included if needed. 

In Paper I we have shown that only the C2 and CN 
(3,0) band of HD 56126 are optically thin. In order to fa- 
cilitate the analysis, we have limited ourself to the C2 (3,0) 
and CN (3,0) bands. For each molecular band we inves- 



tigated whether the different branches for a given lower 
level gave the same column density. Within the errors of 
the determination of the equivalent width this was indeed 
the case. In case the lines would have been optically thick, 
this would have shown up as vertical scatter in the rota- 
tional diagrams (Fig. 5, and 6). This is not observed, and 
it confirms that the lines are optically thin. 

Recently Bakker et al. (1996d) studied the C2 and CN 
bands of HD 56126 by means of a curve of growth analysis. 
They find that the (3,0) bands are indeed close to being 
optically thin. However, in this work any errors which are 
introduced due to the assumption that the lines are opti- 
cally thin are negligible with respect to the errors in the 
measurements of the equivalent width. 

CH+ absorption and emission is assumed to be op- 
tically thin based on the narrow correlation of the data 
points in the rotation diagram (Fig. 7), and the fact that 
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CN (3,0) IRAS 04296 + 3429 
J,,, = 25±6 K 
log N„i,3 = 15.91 ±0.10 p cm ^ 
: Vq=-63.6±0.6 km s"' 


CN (3,0) IRAS 051 13+1347 
= 19+2 K 
log N^t,s = 15.89+0.10 p cm ^ 
: Vq= -3.6+1.1 km s"' 


CN (3,0) HD 56126 
T,,j = 25 + 8 K 
log N„t,3 = 15.51 +0.10 p cm"^ : 
: Vo= 78.3+0.2 km s"' : 


UN [Z^V) IKAb UoUUo — ZOOb 

T,„j = 38 ±2 K 
log N„t,3 = 15.48±0.10 p cm ^ 
: 22.3±0.2 km s"' 
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L T,„, = 30+4 K 

log N„b3 = 1 6.34 + 0.10 p cm ^ 
: Vo=-55.0 + 0.3 km s"' 
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L T,,j = 24+17 K 
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: Vg=-51. 8+0.1 km s"' 


CN M 0"! IRAS 9^^04 + 6147 
L T,„j = 43 + 14 K : 
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Fig. 6. Rotational diagrams of the CN Red System (3,0) band (if not detected the (2,0) band) for the stars showing CN 
absorption. Note that the energies have been computed relative to N" = 0. 
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Fig. 7. Rotational diagrams for the stars showing the CH"*" (0,0) absorption or emission. For CH"*" emission the energy levels 
are relative to A^II J' = 1. 
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Fig. 8. The rotational diagrams of the C2 Phillips (2,0) and (3,0) and CN Red System (2,0) absorption bands in the spectrum 
of IRC +10216. The lower panels gives the part of the spectrum on which we have based our firm identification of C2 and 
tentative identification of CN (see text for details). 



different branches (P, Q, and R) do yield comparable col- 
umn densities for the J" (or J') level they originate from. 
The relative rotational diagram for CH+ emission is rel- 
ative to the energy level of the A^II J' = 1 level. This 
energy level is at 42400 K (over forty thousand degrees !) 
above ground level. 

If the population distribution over the rotational lev- 
els follows a Boltzmann distribution, then the rotational 
diagram should be a linear relation for optically thin lines 
(r < 1). The slope of the curve is inversely proportional 
to the rotational temperature Trot and the offset is given 
by the natural logarithm of N^o\/Qr, where N^noi is the 
column density of the molecule and Qr the partition func- 
tion. If the rotational temperature is not equal to the 
kinetic temperature and different for each pair of levels, 
then the population distribution over the rotational lev- 
els is non-Boltzmann and the rotational diagram will be 
non-linear. From linear least-squares fits to the rotational 
diagrams, we find typical values of T^ot = 43 — 399 K, 
logA/'moi = 14.90 - 15.57 cm-2 for C2 (3,0) (or (2,0)), 



Trot = 155-202 K, logAUoi = 14.35 cm'^ for CH+ (0,0), 
and Trot = 18 to 50 K, logTVmoi = 15.03 to 16.47 cm'^ 
for CN (3,0) (or (2,0)) with an average particle ratio of 
Ar(CN)/Ar(C2) = 4.0±3.0 

Since the rotational temperature is determined using 
all available data points, it is an average temperature. For 
C2, use of only the low J" levels (which are most sensitive 
to collisional (dc-)cxcitation) results in a lower rotational 
temperature, which will be closer to the kinetic tempera- 
ture of the gas. On the other hand, using only the high- J" 
levels (which are sensitive to radiative pumping) results in 
a rotational temperature which will be closer to the color 
temperature of the exciting stellar radiation field. 

The column densities arc computed using the partition 
function, which depends on the average rotational tem- 
perature. In Paper I we investigated the accuracy of this 
method and found that the column densities are slightly 
underestimated for super-thermal C2, and slightly over- 
estimated for sub-thermal CN. This introduces a relative 
error of about 10%. A slightly more accurate determina- 
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tion of the molecular column density could be obtained 

by adding the column densities of all observed J" levels 
and extrapolating to the unobserved levels. The oscillator 
strengths are not very well known and we adopt an ab- 
solute error of 15%. As a consequence, we will adopt an 
absolute error of 15% and a relative error of 10% to the 
molecular column densities. 

The detection of C2 and CN in the same spectrum 
of the same star with significantly different rotational 
temperatures can be attributed to the difference in the 
dipole moment of the two molecules. C2 is a homonu- 
clear molecule (the electronic configuration is spatially 
symmetric with respect to the rotational axis, and the 
nuclear spin of ^^C is 7 = 0). Thus, without a permanent 
dipole moment the selection criteria forbid pure rotational 
and rotational-vibrational transitions in the ground elec- 
tronic state (except through weak quadrupole transitions). 
A homonuclear molecule has therefore no efficient tran- 
sitions available to emit a photon. The higher-J" levels 
will be overpopulated relative to a Boltzmann distribu- 
tion at the local kinetic temperature, because they can 
be populated by radiation with a color temperature larger 
than the local kinetic temperature (optical pumping) . The 
molecule is super-thermally excited. In Fig. 5 this shows 
up as a flattening of the curve for J" > 20 levels. CN is a 
heteronuclear molecule and selection criteria allow dipole 
transitions: energy can be easily released and the molecule 
cools to sub-thermal values if the density is below the den- 
sity for coUisional excitation. 

Remarkably, there seems to be a relation between 
the molecular column density and the expansion veloc- 
ity. Since CN can only be photodissociated by far UV 
photons (A < 1100 A), it is more sensitive to the inter- 
stellar UV radiation field, this leads to a large spread in 
the observed CN column densities. A linear least-squares 
fit (excluding IRAS 08005-2356, IRAS 05341-^0852, and 
IRC 10216) yields: 

logA^ (C2) = 3.53 X IQ-^ X wexp (C2) + 14.98 (3) 
logiV (CN) = 6.19 X IQ-^ X -wexp (CN) + 15.3 (4) 

with correlation coefficients of 0.89 and 0.53, respec- 
tively. Column densities and velocities are in units of cm~^ 
and km s~^ respectively. A detailed discussion of this re- 
lation will be given in § 4, where we will argue that it 
reflects differences in intrinsic carbon abundance. 

3.2.3. Time scales of AGB ejecta 

The far-infrared radiation (IRAS data) is due to thermal 
radiation of dust in the AGB ejecta. From the IRAS 12 
to 25 /um flux ratio, /„ (12/im) / fy (25/im), we derived the 
color temperature of the innermost dust using the method 
described in the IRAS Explanatory Supplement (1986) 
and Table Suppl. VI. C. 6. (Table 1). For stars with a near- 
infrared excess (e.g., a circumbinary disk) this method 



cannot be applied, and we cannot calculate the mass-loss 

rate as described here. 

Trams (1991) has fitted an optically thin dust model 
to a dozen post-AGB stars and found that for all stars 
the dust emissivity efficiency, Q (v) = Qq {v /vqY , has an 
index parameter p =1. Following Sopka et al. (1985), the 
dust inner radius (ro) is given by: 

ro _ / Tdust y(^) 

V0.6Tefi7 ■ ^ ' 

with Tfjust and Teff the Blackbody dust temperature and 
the stellar effective temperature, respectively. Adopting 
logL* = 3.86 log L© for a 0.6 Mq post-AGB star (Blocker 
1995) the stellar radius, can be calculate. Typical dis- 
tances derived are of the order of 10^^ cm (Table 6). Tak- 
ing the dust inner radius and the expansion velocity, we 
can subsequently estimate the time since the star left the 
AGB (Eq. 6, Table 6) and the average annual increase in 
effective temperature (Eq. 7, Table 6): 

ipost-AGB = — — (6) 
^exp 

AT _ Tcff - Tagb 

At tpost-AGB 

while taking Tagb = 3500 K. 

The stars in our sample left the AGB typically 300 
years ago and they have a typical annual effective temper- 
ature increase of 5 K. Since an abundance analysis allows 
the determination of the effective temperature as accurate 
as 100 K (for supergiants), this suggests that we expect 
to see change in the effective temperature and conditions 
of the molecules on a time scale of about 20 years. 

3.2.4. Mass-loss rate 

Here we will attempt to derive the mass-loss rates from 
the column densities and expansion velocities in the AGB 
ejecta. In order to do so, we will make two important 
simplifications: 

(i) Since there is no simple, unique way to determine 
the distance from the star to the molecules, we will as- 
sume that the moleciiles reside where dust is present. In 
the case of the extended envelope of AGB stars and in the 
circumstellar shells of post-AGB stars (F to G-type stars), 
simple molecules are formed out of complex molecules by 
photodissociation by the interstellar radiation field. How- 
ever, for circumstellar shells the abundance of a simple 
molecule will peak at a distance larger than the dust in- 
ner radius. Therefore, the assumption that the molecules 
are formed at the dust inner radius and beyond underes- 
timates the derived mass-loss rate. 

(ii) Our second assumption is that the molecular abun- 
dances are the same for each star, resemble those of the 
AGB star IRC +10216, and do not change with distance 
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Table 5. Derived stellar and dust parameters assuming logL, = 3.86 log Lq. 
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HD 213985 
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27 
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near-ir excess 


IRAS 23304+6147 


5200 
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170 


1440 
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to the star. Prom our finding of a relation between the ob- 
served molecular column density and the expansion veloc- 
ity we know that this is not the case. Furthermore, we take 
the standard abundances to be the computed peak abun- 
dance for IRC +10216. No accurate prediction for molec- 
ular abundance of CH"*" in AGB ejecta was found in the 
literature, while, in the interstellar medium the large ob- 
served abundance of CH+ is still not well understood (see 
Gredel et al. 1993 for a recent summary). 

Using assumptions (i) and (ii) we will derive an equa- 
tion for the mass-loss rate in a rather simple manner. 
Starting with the general formula for conservation of mass, 
M(r) = 4:'jTr^p{r)vexp{r), and assuming a constant mass- 
loss rate over time, p{r) = po (ro/r)^, we find for the den- 
sity at the dust inner radius:: 



Po = 



ro 



mol 



(8) 



where p. is the average hydrogen (molecular, atomic, and 

ionic) particle mass in units of nip and A^moi the molecu- 
lar column density. For the densities and temperatures 
expected in the AGB ejecta most hydrogen will be in 
the form of H2 (/i = 2.0). In the model of Cherchneff 
et al. (1993) the molecular particle abundances relative 
to hh = n(H) + 2 X n(H2) are Xc^ = 4 x IQ-^ and 
^CN = 3 X 10~^. Combining Eqs. 5 & 8 yields the ex- 
pression for the mass-loss rate: 



with all parameters in cgs units except M which is in Mq 
yr-^: 

The mass-loss rates (assuming logi* = 3.86 log Lq) de- 
rived in this manner (with the theoretical abundances of 
IRC 10216) are of the order of -6.2 < logM < -4.1 (Ta- 
ble 5), which is of the same order of magnitude as de- 
rived from the CO emission and the IR excess. In view of 
the important assumptions made in the calculation of the 
mass-loss rate we adopt an estimated error of one order of 
magnitude. We want to stress again that the assumption 
made in calculation the mass-loss rate from the molecular 
absorption lines only allow an order of magnitude esti- 
mate of the mass-loss rate. We want above all to show 
that these molecular lines do allow the determination of 
the mass-loss rate given the right inner and outer radius 
and the molecular abundance. If the mass-loss rate is know 
(from e.g., CO or infrared) the process can be reversed and 
this would yield the molecular abundance. 

A linear least-squares fit (excluding IRAS 08005-2356, 
IRAS 05341+0852, and IRC +10216 for C2 and CN and 
HD 235858 for CN) yields: 



logM (C2) = 0.11 X Wexp (C2) - 6.71 

logM (CN) = 0.13 X wexp (CN) - 5.86 



(10) 
(11) 



with correlation coefficients of 0.95 and 0.92, respectively. 
The mass-loss rate and velocity are in units of [Mq yr~^] 
and km s~^ respectively. 



M = 




3.71 X 10 



-29 -^mol rp 



-(^) 



-^mol 



dust- ''exp- 



(9) 



4. Discussion 

4-.1. General discussion 

We discuss the following points: 

1. Within our sample of stars, eleven stars exhibit C2 
and CN absorption, one CII+ absorption, and one CH+ 
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emission. All post-AGB stars with the Phillips bands in 

absorption also show the Swan bands in absorption. Since 
the a^IIu state is 612 cm~^ above the ground state, cor- 
responding to J" = 18 in the X^S+ state, it is likely that 
this state is populated by optical pumping. A study of the 
Swan band will be subject of a separate paper. 

2. We note that all the detected molecular bands are 
clearly detected, while the non-detections have no molec- 
ular bands or their strength is below our detection limit. 
It therefore seems that there is no star with intermediate 
strength molecular bands. Whether this is due to small 
number statistics or gives clues about circumstellar chem- 
istry is to be evaluated in follow-up research. A sugges- 
tion for the latter option could be that the CSE is rather 
clumpy and that only clumps having a density higher than 
a critical value are able to sustain the presence of these 
molecules. On the other hand one might suggest that these 
molecular bands are rather common, but that our sample 
is too small to include stars with molecular bands of in- 
termediate strength. 

3. All stars exhibiting the unidentified 21 /xm in- 
frared features show the presence of C2 and CN ab- 
sorption, but this is not true for the reverse. The ex- 
ception, IRAS 08005-2356, has the highest column den- 
sity of C2 in our sample, and exhibits OH maser emis- 
sion at the same velocity. At the time of the observation, 
IRAS 05341-1-0852, was the only 21 sources not in- 
cluded in our sample. A recent study on the optical spec- 
trum of this star by Reddy et al. (1997) has indeed show 
the presence of circumstellar C2 and CN absorption. 

We therefore predict that the recently found Post- 
AGB 21 /im sources (IRAS 22574-^6609, Hrivnak & Kwok 
1991a, IRAS 15553-5230 and IRAS 17195-2710, Henning 
et al. 1996, SAO 163075, Justtanont et al. 1996) will show 
C2 and CN absorption. Whether the other new candidate 
and waiting list 21 fim stars (YSO and HII regions) listed 
by Henning et al. show molecular absorption, needs to be 
investigated. But if they do, then there is a one-to-one rela- 
tion between the occurrence of the 21 /xm emission feature 
and the presence of C2 and CN absorption. 

We did not found any clear correlation between the pa- 
rameters listed in Table 5 and the strength of the 21 /xm 
feature (Justtanont et al. 1996, Henning et al. 1996). 
Though there is a suggestion for an anti-correlation be- 
tween the N{CN)/N{C2) and /(20Atm)//(18/xm). Unfor- 
tunately, small number statistics do not allow any firm 
conclusion. 

4. The presence of C2 and CN or CH+ seems to be 
correlated with the dust color temperature: T^ust < 300 K 
and Tdust > 300 K respectively. C2 and CN are correlated 
with cold dust (far-infrared excess) and CH+ with hot 
dust (near-infrared excess). This idea is supported by the 
rotational temperature of the heteronuclear species: Trot ~ 
34 and 200 K for CN and CH+ respectively 

5. The stars in our sample typically left the AGB 
300 years ago and evolved to the left in the HR dia- 



gram with a typical annual temperature increase of 5 K. 
IRAS 08005-2356 has an average annual temperature in- 
crease of 70 K per year. We predict that we can see evo- 
lutionary changes of the star and the CSE in the next 20 
years. 

6. Traditionally, the expansion velocities of the AGB 
ejecta are determined from the full width of the CO(J = 
1^0) millimeter line profile. Since CO is abimdant 
throughout the whole AGB ejecta, the profile contains 
the integrated emission. In the case of expansion veloc- 
ities derived from line absorption, e.g., C2, CN, and CH"*", 
the line-forming region is along the line-of-sight towards 
the stellar photosphere. A second important difference is 
that the molecules observed in absorption are only present 
in a thin shell of material where the conditions are such 
that the molecule has a large abundance. This allows us 
to study the AGB ejecta at different radii, using different 
molecules. Unfortunately, from a theoretical study on the 
extended envelope of IRC-M0216 (Cherchneff et al. 1993), 
is seems that the line-forming regions of C2 and CN are 
almost the same. The upper panels of Fig. 4 clearly show 
the correlation between CO and the C2 and CN expan- 
sion velocity, the relations are given in Eqs. 1 and 2. This 
unambiguously proves that the molecular absorption lines 
originate in the AGB ejecta and that these are thus cir- 
cumstellar. 

7. The lower panels of Fig. 4 show that the C2 column 
density increases as a function of expansion velocity (Eqs. 
3 and 4). To interpret this result, it is important to realize 
that the line-forming region of C2 is only a thin shell. The 
inner radius is determined at a critical dust column density 
such that stellar and/or interstellar photons can penetrate 
and photodissociate complex molecules like C2H2 (via 
C2H) and HCN into C2 and CN, while the outer radius is 
determined in a similar way for the photodissociation of 
C2 and CN to individual atoms by the interstellar and/or 
stellar radiation field. First, we will address the question of 
which radiation field (stellar, interstellar, or circumstellar) 
is responsible for the excitation and the photodissociation 
of C2 and CN. Fig. 9 shows the interstellar radiation field 
(Draine (1978) for A < 2000 A and van Dishoeck & Black 
(1982) for A > 2000 A), and a T^s = 6500 K stellar model 
(Kurucz 1979) with a Tdust = 200 K dust shell. The ra- 
diation field of the star is scaled to a post-AGB star of 
logL = 3.86 log L0 observed at a typical dust inner ra- 
dius of To = 2354 R*. For wavelengths A > 1400 A the 
stellar radiation field dominates, while for A < 1400 A the 
interstellar radiation field dominates. Since C2 is pumped 
by optical radiation (1300 A to 1.1 /xm), the stellar radia- 
tion field is responsible for the excitation of this molecule. 
Since the CN molecule can only be photodissociated by 
photons with A < 1100 A , the interstellar ultraviolet ra- 
diation field photodissociates CN. The C2 molecule is pho- 
todissociated by photons with A < 2000 A and can there- 
fore be dissociated by both the ultraviolet interstellar and 
stellar radiation field. 
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Fig. 9. The interstellar radiation field (solid straight line) and the stellar radiation field at To ~ 2719 R* of a Tbs = 6500 K, 
log 3 = 1 Kurucz model. A typical infraxed excess has been visualized by a Tdust = 200 K black-body (scaled arbitrarily) . Optical 
pumping of C2 uses photons between 1300 A and 1.1 /xm (horizontal dashed line), while photodissociation of C2 and CN needs 
photons with A < 2000 and 1100 A, respectively. 



We have clearly found a relation between the expan- 
sion velocity and the column density of C2 and to a lesser 
extend of CN. The sign of the slope is inconsistent with 
models which assume a constant carbon abundance and 
mass-loss rate: these models predict a decrease in column 
density for increasing expansion velocities. 
Alternatively we propose two models: 
Our favorite model (model I): The radiation pressure effi- 
ciency {Qrp see Tielens 1983) scales with the carbon abun- 
dance. Since the momentum equation gives a positive cor- 
relation between Q^.^ and the acceleration dvg/dr, this 
would account for the spread in the observed expansion 
velocities. At the same time the spread in carbon abun- 
dance in the ejecta will be observable as a spread in the 
observed column densities in a way consistent with our 
findings. We ran some simple models and found that only 
a small increase of Q^p gave rise to a significant change 
in expansion velocity. Additional and much more detailed 
modeling is needed to determine the exact relation be- 



tween Uoxp, log A'', and Q^p. 

We propose an alternative model (model II): The higher 
the initial mass of the star, the higher the mass-loss rate 
and terminal velocity of the wind when the star termi- 
nates the AGB evolution (Barnbaum et al. 1991). Since 
an increase in mass-loss rate increases the column den- 
sity, but an increase of the expansion velocity decreases 
it, the model can only explain the observation if the in- 
crease of mass-loss rate dominates over the increase of the 
expansion velocity. Since the number of thermal pulses 
increases with the mass of the progenitor, we also expect 
an increase of carbon abundance which would work in fa- 
vor of this model. Therefore model II included the effect 
discussed for model I. 

Since the model with the least free parameters is 
stronger than a more complex model, we favor model I 
and argue that a high C2 (and CN) abundance reflects a 
high carbon abundance (and the C/0 ratio), which in turn 
implies that the stars have experienced different amounts 
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(number or efficiency) of third dredge-ups. This hypoth- 
esis implies that there shoiild be a relation between the 
carbon abundance of the star and the observed column 
density of circumstellar C2. 

We further note that there seems to be a flattening 
(possibly due to saturation of the absorption lines, flat 
part of curve of growth) in the molecular column densities 
for Wexp ^ 30 km s"'^ , but whether this is real or an artifact 
of small number statistics is not clear at this moment. 

8. There should be a reasonable explanation why the 
relation between log N and Ucxp for CN is not as well de- 
fined as for C2. Although the CN molecular spectra are 
more complex than the C2 spectra, and a unique identi- 
fication of an absorption line with a single transition is 
not always possible, it is unlikely that this can account 
for the large scatter around the mean relation. Noting 
that all stars with m„ > 12 have N{CN)/N{C2)> 3, 
while all stars with to„ < 12 have N{CN)/N{C2)< 3, 
we suggest that the CN abundance might be affected by 
the distance of the star above the Galactic plane, and 
the local interstellar ultraviolet radiation field. It might 
also be that the fainter objects have a larger circumstellar 
reddening and therefore the molecules are more shielded 
against the stellar/interstellar radiation field. This would 
increase A^(CN). However, information about the circum- 
stellar reddening of these objects is not available. 

9. Mass- loss rates have been derived from the observed 
column densities by making two very important assump- 
tions. The basic conclusion is that the mass-loss rates de- 
termined are of the same order as those obtained from CO 
emission and infrared measurements, but do not allow a 
detailed comparison between different objects. 

10. Fedcrman et al. (1994) presented a relation be- 
tween the C2 and CN column density observed in dif- 
fuse interstellar clouds. We have investigated if this rela- 
tion could be extrapolated to the conditions prevailing in 
the CSE of post-AGB stars (Fig. 10). We find (excluding 
IRC -1-10216, and IRAS 05341-1-0852): 

logAr(CN) = 1.153 X logiV(C2) - 1.863 (12) 

with N in units of cm~^ with a low correlation coefficient 
of 0.43. Although both in the interstellar and circumstel- 
lar case C2 and CN are primarily destroyed by photodis- 
sociation by the interstellar radiation field, the formation 
chemistry of these species is assumed to be very different 
in interstellar clouds. 

11. We have searched for the lower v' bands (since 
they are the strongest); we did not detect a single iso- 
topic line of C2, CN, or CH"*". Depending on the spectrum 
the detection limit is about 5 to 10 mA, which gives a 
typical lower limit of ^^C/-^^C > 20. ^^C is produced by 
the CN-cyclc and not by the triplc-a process. A high iso- 
tope ratio indicates that the carbon enhancement is due 
to the third-dredge up (convection reaches the deeper He- 
burning shell). This is rather surprising since one would 



naively expect that the material from the outer shell (H- 

burning) can be dredgcd-up easily. The high isotope ra- 
tios of these stars is consistent with carbon stars (AGB 
stars) to be the progenitor of the stars studied (^^C/^^C 
« 40 — 80, Lambert et al. 1986). The third dredge-up only 
occurs on the thermal pulsating AGB phase for stars with 
.^MS < 5 Mq {Z = 0.001): this sets an upper limit on 
the star's initial main sequence mass. For higher initial 
masses, Hot Bottom Burning (HBB) will prevent the for- 
mation of a carbon rich post-AGB star (Boothroyd et al. 
1993). 

12. Ultraviolet: Lambert et al. (1995) have obtained 
HST spectra of diffuse interstellar clouds, showing the 
Mulhkcn (0-0) (2313 A) and F-X (0-0) (1342 A) bands of 
C2 in absorption. CO (e.g., 1509 and 1478 A) and excited 
H2 (e.g., 1108 and 1092 A) are also prominent absorbers in 
the UV, since the column densities in the AGB ejecta are 
about two orders of magnitude larger than for interstellar 
clouds, the ultraviolet spectra of post-AGB stars as mea- 
sured, e.g. with the HST and lUE, should be dominated 
by molecular (absorption) bands. 

13. Visual: We predict the presence of the CN Violet 
system and the C3 Swing system. Many molecules have 
been observed in the optical spectra of Comets. Though 
these are 0-rich environments we would expect some of 
these molecules also to be present in C-rich environments. 
We therefore suggest the presence of the Merril-Sandford 
SiC2 bands (4640 and 4977 A), CH (3130-3150 A) and NH 
(3358 A). 

14. Infrared: Molecules have many transitions in the 
infrared. For the stars in our sample we might expect 
molecular absorption and emission in the infrared of e.g., 
H2 (2-6 Mm), HCN (3.4 ^m), HCN (14.1 /xm), C2H (27.1 
fim), C2H2 (13.6-13.8 iim). 

13. Millimeter: Only for the strongest transitions of 
abundant molecules are lines in the optical and ultraviolet 
spectra observed. Many molecular transitions have been 
observed in the (sub-)millimeter and radio for IRC +10216 
(for details see Kawaguchi et al. 1995), but very little work 
has been done to look at molecules in these wavelengths 
for post-AGB stars. However, to determine accurate abun- 
dances and to understand the circumstellar chemistry of 
detached dust shells, observations such as those presented 
here are of crucial importance. 

4-2. Discussion of individual objects 

IRAS 04296-1-3429: Very typical in its circumstellar CN 

and C2 absorption. 

IRAS 05113-1-1347: Very typical in its circumstellar CN 
and C2 absorption. Our spectra show photospheric CN 

bands. 

IRAS 05341-1-0852: This star has an overabundance of 
s-process elements and carbon (Reddy et al. 1997) which 
clearly suggests that this is a post-AGB stars. Because of 




Fig. 10. A logarithmic plot of N(CN) versus N(C2). The thick dashed line represents the interstellar data points from Federman 
et al. (1994), while the circumstellar data points are determined in this study. Is seems that the circumstellar points follow the 
interstellar trend. 



the noisy spectra the data on C2 and CN is not of the 
same quality as those for the other stars. 
HD 44179 (The Red Rectangle): We confirm the pres- 
ence of CH+ (0,0) and (1,0) in emission at a projected ex- 
pansion velocity of 3.4 ± 2.0 km s^^ with Tj-ot = 202 ± 16 
K, and the non-detection of C2 and CN. CH^ in the Red 
Rectangle was first detected by Waclkens et al. (1992) 
and later identified by Balm & Jura (1992) and Hall et 
al. (1992) and is so far the only object which shows CH+ 
in emission. The latter authors estimated the rotational 
temperature to be Trot = 120 ±50 K. We prefer our deter- 
mination of the rotational temperature since many more 
linos are included and the error on the equivalent width is 
significantly smaller. 

Although we did not detect CN nor C2, the latter 
molecule is observed in emission in the reflection lobes 
of the Red Rectangle (Sarre 1996). Schmidt et al. (1980) 
found numerous narrow features in the emission spectrum 
of the lobes, but were not able to make a proper identi- 
fication. Glinsky et al. (1996a) review these unidentified 



molecular bands and argue that they are due to phos- 
phorescence of C3. Glinsky et al. (199b) report on their 
detection of the spin-forbidden Cameron bands of CO in 
the ultraviolet. 

The absence of C2 and CN absorption could be related 
to the low CO/I(60 /xm) ratio (van der Veen et al. 1993). 
Since CO is underabundant and the star itself is extremely 
metal-depleted this might be the result of the same pro- 
cess: condensation of gas and molecules on circumstellar 
dust grains. Jura et al. (1995) detected weak CO millime- 
ter radiation at v^,_q = 18.9 ± 2.0 km s"-'^ (being the he- 
liocentric stellar velocity derived from CO emission) with 
Vexp ~ 6 km s~^. This is close to the CH+ central velocity 
and suggests that CH"*" emission is exactly on the system 
velocity. In the spectra discussed CH+ emission lines are 
not resolved. New observations (Bakker et al. 1996d) at 
R « 120, 000 resolve the line profiles of the stronger lines 
with a FWHM ^ 8.5 km s'^ and FWFM « 15 km s~i 
(full width at which the emission falls below the detection 
limit). Since the formation of CH+ is probably due to 



Eric J. Bakker et al.: Circumstellar C2, CN, and CH"'" in the optical spectra of post-AGB stars 



19 



shocks in a circumbinary disk (the emission comes from 

levels 42000 K above the ground level), the FWFM is 
an indication for the turbulent and Kepplerian velocity of 
the line forming region. The interstellar CH+ abundance 
is not very well understood. Early models (e.g., Elitzur & 
Watson 1980) have suggested that CH"*" might be formed 
efficiently in shocked regions. Although such models are 
currently not favored for the CH"*" production in interstel- 
lar clouds, it could take place in the complex environment 
of the Red Rectangle. The fact that no absorption com- 
ponent is observed suggests that the line-forming region is 
not very extended and is within the "slit" of the telescope 
(< 2"). 

It is interesting to note that the diffuse interstellar 
bands are in emission (Sarre 1991), although we did not 
find these bands in emission in our spectrum. Roddier et 
al. (1995) have spatially resolved the near-infrared emis- 
sion into two separate emission peaks with an angular 
separation of 0.14". Waelkens et al. (1996) have clearly 
demonstrated that it is a single- lined spectroscopic binary 
with an orbital period of 318 ± 3 days. The latter authors 
argue that the star is not directly observed but only via 
scattered light on the reflection lobes. 
HD 52961: Like the Red Rectangle and HD 213985, the 
photospheric abundance pattern shows evidence for selec- 
tive accretion of circumstellar gas, but unlike these two 
other objects, no CH"*" has been detected. 
HD 56126: The C2 and CN bands have been studied 
in detail in Paper I. It was shown that optical depth ef- 
fects are important and that the excitation of C2 is a 
balance between radiative and collisional (de-)excitation. 
HD 56126 has 21 (Kwok et al. 1989) and 30 /im features 
(Omont et al. 1995) and shows many of the s-process el- 
ements in its spectrum (Klochkova 1995). Oudmaijer & 
Bakker (1994) suggested a pulsation period between 30 
and 96 days. More recently Lebre et al. (1996) discuss 
the possible RV Tauri nature of HD 56126, and they sug- 
gest the presence of four periods (32.9, 27.3, 11.9, and 7.1 
days). We did a similar analysis with the data available 
from Lebre et al. and some of our own data and con- 
firm the presence of four significant minima in the Phase 
Dispersion diagram on the radial velocities of HD 56126. 
However, it seems possible that there is only one period 
of Ppuis = 12.1 days and that the other minima are (sub-) 
harmonics. 

IRAS 08005-2356: Slijkhuis et al. (1991) have shown 

that in the optical spectrum all absorption lines are blue- 
shifted up to velocities of 50 km s~^. They interpreted this 
as a wind spectrum resulting from a huge post-AGB mass- 
loss rate of about 10^^ Mg yr"-*^ with a terminal velocity 
of about 400 km s"^. The Ha profile shows significant 
changes within a time interval of less than two months 
which suggests a strong asymmetric mass-loss. The ACB 
ejecta emit OH maser emission with a maximum (pro- 
jected) velocity of 50 km s"^ (Te Lintel Hekkert 1991), 
and the infrared energy distribution shows the presence of 



a cold dust shell {R > 3100 R*, Tdust < 150 K) and a hot 
dust shell (or disk) at i? > 34 R* with Tdust < 1200 K. 
The overall observational characteristics of IRAS 08005- 
2356 are very similar to those of HD 101584 and suggest 
that the model by HD 101584 of Bakker et al. (1996b), 
consisting of a close binary system with an eccentric orbit 
also applies to IRAS 08005-2356. The secondary triggers 
the high- mass loss of the primary and the stellar wind acts 
as a curtain in front of a star much hotter than F5I. This 
suggests that the optical spectrum probably does not con- 
tain information about the star, but only about the wind. 

The presence of C2 and C3 was also noted by Hriv- 
nak (1995). Using our high-resolution spectra we are able 
to determine the expansion velocity of 50 km s~^. C2 
absorption has the same radial velocity as the blue OH 
maser emission peak. The fact that the OH maser emis- 
sion and C2 are observed at the same velocity is unex- 
pected since the first is a tracer of oxygen rich material 
(H2O dissociation product), while the latter is an indica- 
tor of carbon- rich material (C2H2 dissociation product). 
This suggests that the star has only recently changed from 
oxygen to carbon-rich without changing its terminal wind 
velocity, or that there are separate carbon and oxygen 
rich mass-losing stars in the system, or that the circum- 
stellar chemistry is very unusual. Bakker et al. (1996d) 
have observed the molecular absorption lines at a reso- 
lution of i? w 120, 000 and resolved the molecular fea- 
tures in two separate component with a line splitting of 
5.7±2.0 km s~^. The most likely explanation is that there 
are two photodissociation fronts at different velocities: one 
due to the stellar and one due to the interstellar radiation 
field. We further note that there are about twenty chro- 
mospheric narrow emission lines from neutral and singly 
ionized metals present in our spectrum, which are double 
peaked, consistent with an accretion disk. The presence 
of an accretion disk, and the accretion of material can 
generate a high velocity bipolar outflow. 

We not that this is the only star in our sample which 
shows C2 and CN but not the 21 i^m feature. Very cau- 
tiously we propose that this object might have a 21 /xm 
feature but that it was excluded from detection so far. 
IRC +10216: Much of our knowledge about carbon- 
rich circumstellar environments is based on observing and 
modeling the carbon- rich ACB star IRC +10216. Over the 
years forty different molecular species have been identified 
in the radio and (sub) millimeter spectrum of this star (Lu- 
cas 1992). The mass-loss rate and expansion velocity have 
been determined from CO observations by Huggins et al. 
(1988) at 4+1 X 10"^ Mq yr'^ and 14+1 km s'^ Circum- 
stellar C2 and visible lines of CN have not been reported 
and here we present the first detection of C2 and possibly 
CN at an expansion velocity of 13.8 + 2.0 km s~^ (Fig. 8). 
The spectrum is severely blended and an accurate deter- 
mination of the continuum level and equivalent widths is 
impossible. Taking the local pseudo continuum as the con- 
tinuum level we have determined equivalent widths and 
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derived T^t = 43 ± 6 K, log AT = 14.90 ± 0.10 cm-^, and 
vq = -33.3 ±0.7 km s-\and Trot = 26 ± 11 K, logA^ = 
14.90±0.10 cm-2, and vq = -32.3±0.7 km s^^ for the C2 
(2,0) and (3,0) bands, respectively. The theoretical work 
by Chcrchneff et al. (1993) predicts Xq^ = 4 x 10"^ for a 
distance between 50 and 200 x 10^^ cm. With a mass-loss 
rate of 3 x 10~^ Mq yr~^ and an expansion velocity of 14 
km s^^ (Huggins ct al. 1988) this gives a column density of 
log A^predicted = 15.30 cm~^. Our observations are consis- 
tent with the predictions if we assume that the real contin- 
uum is a factor of three higher than the observed pseudo 
continuum. Alternatively, the C2 abundance of Cherchneff 
et al. (1993) is a factor three too high. 
HR 4049: This truly remarkable object is the most metal- 
depleted star known ([Fe/H]~ -4.8, Waelkens et al. 1991). 
Bakker et al. (1996a) have studied the complete optical 
spectrum at high-resolution and high signal-to-noise ratio 
without detecting a single iron peak line. A strong near- 
infrared excess has been detected (Lamers et al. 1986) 
which is probably due to the presence of a circumbinary 
disk. However, unlike the Red Rectangle and HD 213985, 
no CH+ absorption or emission has been detected in 
HR 4049, which suggests that the post-AGB mass-loss, 
or mass-transfer in the binary system, is less violent and 
does not result in the presence of a shocked region of cir- 
cumstellar gas. 

HD 161796: An 0-rich supergiant believed to be a post- 
AGB star. No molecules have been detected in the optical 

spectrum of this star. 

IRAS 20000+3239: Not much is known about this ob- 
ject except that it exhibits the unidentified 21 /xm (Kwok 

et al. 1995) and 30 /im (Omont ct al. 1995) features. In 
the sample of stars with C2 and CN absorption lines this is 
the star with the highest N{CN)/N{C2) ratio of 11.2, the 
hottest CN rotational temperature of Tj-ot = 50 ± 12 K, 
while it has a typical expansion velocity of 12.8 km s"^. 
The CN column density relative to C2 is about twice as 
high as the average. Our spectra show photosphcric CN 
bands. We also note that the CO absorption in the K-band 
is the strongest in the sample of stars observed (Hrivnak 
et al. 1994). 

AFGL 2688 (The Egg Nebula): This object is ob- 
served as a highly reddened central star surrounded by 

a torus (possibly a close binary system) with two bright 
lobes (reflection nebulae) at the equatorial poles of the sys- 
tem. Three different stellar winds are observed in the CO 
millimeter line emission (Young et al. 1992): high (HVW), 
medium (MVW) and low velocity wind (LVW) with ex- 
pansion velocities of 100 ± 10, 45 and 22.8 km s~^, re- 
spectively. The HVW is bipolar with a dc-projcctcd wind 
velocity of 360 km 8"^ While the HVW and MVW are 
post-AGB winds, the LVW is the stellar wind when the 
star was on the AGB (now the AGB ejecta). From model- 
ing of the CO millimeter emission Young et al. found that 
the AGB ejecta contain 0.7 Mq fed by an AGB mass-loss 
rate of M = 1.5 x 10"^ M© yr'^ 



C2 emission was first reported by Crampton et al. 

(1975) and later confirmed by Cohen & Kuhi (1977). From 
optical spectropolarimetry, Cohen & Kuhi have shown 
that C2 is in absorption in the reflected, polarized, light 
of the lobes and has an unpolarized emission component, 
which is attributed to emission within the "slit" . Cramp- 
ton et al. also report C3 absorption. 

We have observed the brightest lobe (north lobe) of the 
Cygnus Egg Nebula and found C2 in absorption with an 
expansion velocity of 17.3±2.0 km s~^. This indicates that 
C2 (and CN) is formed in the LVW. However, we did not 
observe any emission, which is probably due to the smaller 
"slit length" used, so that only reflected light from the lobe 
is observed and no emission from the surrounding gas. 
Although reported by Cohen & Kuhi (1977), our spectra 
do not show the presence of the SiC2 Merrill-Sanford band 
(4977 A) in absorption. 

We did not detect the ^^CN Red System (1,0) band 
which places a lower limit of ^^C/^^C > 19. This is con- 
sistent with the the isotope ratio of ^^C/^'^C « 20 found by 
Wannier & Sahai (1987) in the slow wind. The fast wind 
has an isotope ratio of ^^C/^^C w 5 (Jaminet & Danchi 
1992). Jaminet & Danchi foimd a self absorption feature 
in the CN millimeter lines at the same velocity we find 
the optical absorption lines. This opens the possibility to 
study the CN molecule by both its optical absorption lines 
and the millimeter emission lines. 

IRAS 22223+4327: Together with HD 235858 these are 

the only two sources in our sample which show the CO 
first overtone in emission. This is attributed to coUisional 
excitation of the CO molecule in the circumstellar envi- 
ronment. C2 and C3 detections are reported by Hrivnak 
(1995). We confirm the presence of C2 absorption and add 
CN to the list. As a result of the large number of observed 
molecular lines of C2 the rotational temperature and col- 
umn density are among the highest in our sample. 
HD 235858 (IRAS 22272+5435): This photometric 
variable has an impressively strong far-infrared excess due 
to a detached dust shell. Together with the presence of 
strong absorption lines from s-processed elements this is 
clearly a post-AGB star (Zacs et al. 1995). C2 and C3 ab- 
sorption in the optical spectrum were first noted by Hriv- 
nak & Kwok (1991b), while here we add CN to the list 
of detections. Very strong photosphcric CN Red System 
bands lines (Bakker et al. 1996d) and narrow line circum- 
stellar CN Red System bands are observed simultaneously 
(e.g, 7895 to 7905 A). Hrivnak et al. (1994) reported on 
the change from CO first overtone absorption to emission 
in less than two months. Since CO first overtone bands 
are formed in hot gas (Tgas ^ 1000 K) it is not expected 
that circumstellar C2 or CN will follow the variations of 
CO. HD 235858 differs from the other objects showing C2 
and CN absorption in that the CN abundance is the low- 
est in our sample. Since we have argued that the large 
spread in observed column densities of CN is probably 
due to differences in the UV radiation field it seems likely 
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that HD 235858 is exposed to a stronger interstellar UV 

radiation field than average. 

HD 213985: The overall energy distribution is very much 
like that of HR 4049. Differences occur in the ultraviolet 

because of different circumstellar extinction laws, and in 
the far-infrared due to the presence of cool AGB ejecta. 
The infrared energy distribution of HD 213985 can be 
modeled using two Blackbodies: one component at a dis- 
tance of 22 iZ* with Tdust = 1250 K, probably a cir- 
cumbinary disk, and a second component at 220 with 
Tdust = 350 K, the AGB ejecta (Waelkens et al. 1987). 
Here we confirm the detection of CH"'" in absorption at 
an expansion velocity of 6.7 ± 2.0 km (Waelkens et al. 
1995). If CH"*" is formed in the AGB ejecta at a distance 
of 220 (with i?, = 50 Rq) then HD 213985 left the 
AGB only 43 years ago. This would mean that the effec- 
tive temperature has increased from 3500 K to 8500 K in 
43 years giving an average annual increase of 100 K. Since 
this is not observed, CH"*" cannot be formed in the AGB 
ejecta but is probably formed much closer to the star: a 
shocked region in a slowly expanding circumbinary disk. 

We note that in our spectra of HD 213985 there is 
evidence for line splitting very similar to those observed for 
W Virginis stars (e.g., ST Pup and V29, Gonzalez 1993). 
Most noticeable the Till Unes at 4549.61 and 4563.77 show 
not only a narrow absorption feature on the red wing, but 
also the Fell displayed in Fig. 3. 

BD +39°4926: This star is almost a twin of the central 

star of the Red Rectangle, with the difference that no 
infrared excess has been detected for this object. Unlike 
the Red Rectangle and HD 213985, no CH""" absorption 
or emission has been detected. This combined with the 
absence of an infrared excess suggests that the post-AGB 
mass-loss, or mass-transfer in the binary system, is less 
violent and does not result in the presence of a shocked 
region of circumstellar gas. 

IRAS 23304+6147: Very typical in its circumstellar CN 
and C2 absorption. 

5. Conclusions 

We have explored a new technique to study the physi- 
cal and chemical conditions of the AGB ejecta by looking 
at molecular absorption and emission lines in the optical 
spectra of post-AGB stars. We find that all stars exhibit- 
ing the unidentified 21 fim feature have C2 and CN ab- 
sorption. Stars which show C2 and CN do not show CH"*" 
absorption. The presence of C2 and CN is correlated with 
the presence of cold dust (Tjust < 300 K), while CH""" is 
correlated with the presence of hot dust (Tdust > 300 K). 

The expansion velocities determined from the molecu- 
lar absorption lines are in very good agreement with the 
expansion velocities derived from CO millimeter line emis- 
sion. This proves that the molecular absorption lines are 
formed in the AGB ejecta. The absolute heliocentric ve- 
locity of the AGB ejecta in the line-of-sight can be de- 



termined very accurately thanks to the large number of 

molecular lines available in the optical. We find a typical 
error of ~ 0.3 km s""^ in our observations. 

From the observed equivalent widths (see App. A, 
only at CDS) of the absorption lines the rotational tem- 
peratures and column densities can be determined. We 
find that the rotational temperature of C2 is significantly 
higher than that of CN. C2 is super-thermally excited, 
whereas CN is sub-thermally excited. This is consistent 
with the fact that C2 is a homonuclear molecule (with 
7 = 0) and CN a hetcromiclcar molecule, while the pri- 
mary excitation mechanism of C2 is optical pumping by 
the stellar radiation field. A more detailed analysis of the 
excitation of these species can lead to better constraints 
on the physical parameters in the AGB ejecta and to an 
independent determination of the mass-loss rate. These 
points will be further investigated in a subsequent paper 
(Paper HI in preparation). 

Interestingly, we found that the molecular column den- 
sities increase with expansion velocity. This is interpreted 
as due to the fact that carbon-rich dust is accelerated 
to higher velocities by the stellar radiation field. The ob- 
served column densities arc an indicator of the molecular 
abundance. Mass-loss rates are computed which are of the 
same order of magnitude as those found from the IR excess 
and from CO emission lines. In view of the important as- 
sumption made to be able to compute the mass-loss rate, 
we stress that these rates should be cited cautiously. 
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A. Appendix: Equivalent width of C2, CN, and CH+ lines 

MS 5213 

Circumstellar C2, CN, and CH"*" in the optical spectra of post-AGB stars 

Eric J. Bakker, Ewine F. van Dishoeck, L.B.F.M. Waters, and Ton Schoenmaker 

A&A main journal, 1997, volume, first page 

molecular processes — circumstellar matter - stars: AGB and post-AGB — line: identification 

Tables A.f/2/3/4/5 are available in electr onic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (f 30.79. f 28.5) or via 



http://cdswob.u-strasbg.fr/Abstract.html, or from the authors 



Description of parameters: 

B branch identification 

J" the rotational quantum number, 

total angular momentum including spin (Herzberg 1950) 

N" the rotational quantum number, 

total angular momentum excluding spin (Herzberg 1950) 
for CN: J" = iV" - 1/2 (Fa), J" = N" + 1/2 (Fi) 

A Laboratory wavelength of transition in air in A=le-10m 

f{J' J") oscillator strength 

EW equivalent width in mA=le-13m, 

positive values are absorption lines and 

negative values (only for HD 44179) are emission lines. 

Branch identification for C2 A^H^ - XiE+ and CH+ A^Hu - XiS+: 

B = -1 P Branch (AJ = -1 = J' - J") 

B= Q Branch (AJ = = J' - J") 

B= 1 R Branch (AJ = 1 = J' - J") 

Branch identification for CN A^D - X^E+ (after J0r genson & Larsson 1990): 

B= 1 Ri 

B= 2 Qi 

B= 3 Pi 

B= 4 "^Riz 

B= 5 PQ12 

B= 6 °Pi2 

B= 7 R2 

B= 8 Q2 

B= 9 P2 

B= 10 ^R2i 

B= 11 "-Qzi 

B= 12 '3P2i 

Code used to identify a star: 

*042 IRAS 04296-H3429 

*051 IRAS 05113-^1347 

*053 IRAS 05341-^0852 

*441 HD 44179 

*561 HD 56126 

*080 IRAS 08005-2356 

*102 IRC -hl0216 

*200 IRAS 20000+3239 

*268 AFGL 2688 

*222 IRAS 22223+4327 

*235 HD 235858 

*233 IRAS 23304+6147 

*213 HD 213985 

Tables' 

Table'A.l. C2 A^n^ - Phillips (2,0) band. 

Table A.2. C2 A^H^ - Phillips (3,0) band. 

Table A.3.a CN A^H - X^T,+ Red System (2,0) band. 

Table A.3.b Continued: CN A^H - X^E+ Red System (2,0) band. 

Table A.4.a CN A^H - X^T,+ Red System (3,0) band. 

Table A.4.b CN A^H - Red System (3,0) band. 

Table A.5. Continued: CH+ A^H - X^T.+ (0,0) band. 
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Table A.l. C2 A^Hu - X^S+ Phillips (2,0) band. 
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A (air) [AJ 


f(J J ) 


hW [mAJ 
*102 
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8750.850 


A JO T > A 

4.43E-04 


44.3 
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8751.490 


A C\ A T7^ f\ A 

4.24E-04 
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4.80E-04 


117.5 
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10 
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/I 1 1 XT' t~\ A 
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5.76E-04 


77.1 
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"1 A A T7^ 
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-1 
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7.19E-04 




-1 


4 


8773.430 


2.40E-04 


29.6 


1 


18 


8776.611 


3.88E-04 







10 


^TO -I A A 

8780.144 


7.18E-04 




-1 


6 


8782.311 


2.76E-04 




1 


20 


8786.050 


3.85E-04 







12 


8788.561 


7.17E-04 




-1 


8 


8792.652 


2.96E-04 




1 


22 


8796.979 


3.82E-04 







14 


8798.462 


7.17E-04 




-1 


10 


8804.502 


3.07E-04 




1 


24 


8809.410 


3.80E-04 







16 


8809.844 


7.16E-04 




-1 


12 


00-1 ^ r>o/^ 

8817.830 


3.15E-04 







18 


8822.728 


7.15E-04 




1 


26 


8823.379 


3.78E-04 




-1 


14 


8832.682 


3.20E-04 







20 


8837.122 


7.14E-04 




1 


28 


0000 A 

8838.846 


3.76E-04 




-1 


16 


8849.075 


3.24E-04 







22 


K f\ A A 

8853.044 


7.12E-04 




1 


30 


8855.876 


3.73E-04 




-1 


18 


8866.997 


3.27E-04 







24 


8870.519 


7.11E-04 




1 


32 


8874.466 


3.72E-04 




-1 


20 


8886.487 


3.29E-04 







26 


r' r' 

8889.535 


7.10E-04 




1 


34 


8894.623 


3.70E-04 




-1 


22 


8907.545 


3.30E-04 







28 


8910.135 


7.08E-04 




1 


36 


cir\-\ r' A f\f\ 

8916.400 


3.68E-04 




-1 


24 


8930.172 


3.31E-04 







30 


8932.333 


7.06E-04 




-1 


26 


8954.440 


3.32E-04 







32 


8956.139 


7.04E-04 




-1 


28 


8980.336 


3.33E-04 







34 


8981.586 


7.02E-04 




-1 


30 


9007.897 


3.33E-04 







36 


9008.705 


7.00E-04 




-1 


32 


9037.150 


3.33E-04 




-1 


34 


9068.086 


3.33E-04 




-1 


36 


9100.823 


3.32E-04 
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Table A.2. C2 A^Hu - X^S+ Phillips (3,0) band. 



B T A (air) [A] fjJ'J") EW [mlj 











*042 


*051 


*053 


*561 


*080 


*102 


*200 


*268 


*222 


*235 


*233 


1 


6 


7714.578 


2.006-04 


34.9 


60.5 


16.8 


28.0 


85.0 




33.4 


43.3 


52.1 


40.0 


70.7 


1 


4 


7714.947 


2.23E-04 


47.9 


34.0 


27.4 


30.7 


95.0 




33.6 


76.6 


61.9 


17.0 


73.2 


1 


8 


7715.418 


1.96E-04 


22.1 


28.0 


11.6 


24.0 


40.3 




21.9 


25.7 


47.4 


12.0 


51.5 


1 


2 


7716.531 


2.67E-04 


53.9 


23.8 




29.5 


100.0 




57.2 


133.9 


66.9 


36.7 


62.2 


1 


10 


7717.473 


1.91E-04 








14.2 


31.9 




25.1 




29.1 




35.0 


1 





7719.331 


6.67E-04 


58.3 


27.3 


34.2 


30.5 


73.4 


72.0 


55.2 


152.7 


58.9 


59.5 


56.4 


1 


12 


7720.751 


1.87E-04 




13.4 




23.1 


24.7 




25.1 




21.6 




23.3 





2 


7722.096 


3.34E-04 


56.2 


22.1 


44.9 


34.7 


91.7 


44.0 


67.6 


152.4 


75.2 


47.5 


78.2 





4 


7724.222 


3.34E-04 


61.6 


26.1 


30.4 




94.0 


38.0 


63.6 


98.3 


82.9 


59.9 


77.2 


1 


14 


7725.243 


1.84E-04 










12.8 




19.1 




18.0 




12.4 


-1 


2 


7725.822 


6.65E-05 


30.9 


12.4 


9.6 


16.4 


55.8 




33.8 


70.1 


28.9 


22.0 


38.7 





6 


7727.560 


3.34E-04 


41.5 




49.0 




84.8 




25.8 




78.3 


20.0 


71.7 


1 


16 


7730.966 


1.82E-04 




4.3 




11.4 






9.7 




14.0 




8.3 


-1 


4 


7731.666 


l.llE-04 


14.0 


9.1 


26.2 


29.2 


70.1 




39.5 


44.9 


39.2 


45.0 


50.6 





8 


7732.120 


3.33E-04 


22.7 


15.5 


18.2 


29.6 


61.3 




40.5 




64.9 


33.0 


82.0 





10 


7737.908 


3.33E-04 


18.9 


18.3 






44.7 




32.2 




55.0 


65.6 


53.4 


1 


18 


7737.908 


1.80E-04 
























-1 


6 


7738.740 


1.28E-04 


10.3 


11.9 




20.6 


51.8 




32.6 




36.1 


33.5 


51.5 





12 


7744.903 


3.33E-04 










27.3 














1 


20 


7746.098 


1.78E-04 
























-1 


8 


7747.040 


1.37E-04 










47.5 

















14 


7753.144 


3.32E-04 










16.1 














1 


22 


7755.532 


1.77E-04 
























-1 


10 


7756.585 


1.42E-04 










16.0 

















16 


7762.626 


3.31E-04 










11.0 














1 


24 


7766.226 


1.76E-04 
























-1 


12 


7767.372 


1.46E-04 








13.7 


12.9 

















18 


7773.358 


3.31E-04 


1.0 












14.3 




12.0 




7.3 


1 


26 


7778.178 


1.75E-04 
























-1 


14 


7779.431 


1.48E-04 








14.0 










13.7 




12.8 





20 


7785.350 


3.31E-04 


















17.2 




15.7 


1 


28 


7791.405 


1.74E-04 
























-1 


16 


7792.741 


1.50E-04 


















12.0 




14.9 





22 


7798.613 


3.30E-04 


















10.1 




11.2 


1 


30 


7805.908 


1.73E-04 
























-1 


18 


7807.322 


1.52E-04 


















9.8 




9.6 





24 


7813.138 


3.30E-04 


















15.6 




4.5 


1 


32 


7821.720 


1.72E-04 
























-1 


20 


7823.201 


1.53E-04 


















4.7 









26 


7828.979 


3.29E-04 


















8.2 






1 


34 


7838.850 


1.71E-04 
























-1 


22 


7840.387 


1.53E-04 



























28 


7846.122 


3.28E-04 


















10.2 






1 


36 


7857.314 


1.70E-04 
























-1 


24 


7858.889 


1.53E-04 



























30 


7864.595 


3.27E-04 
























-1 


26 


7878.721 


1.54E-04 



























32 


7884.394 


3.27E-04 
























-1 


28 


7899.907 


1.54E-04 



























34 


7905.561 


3.26E-04 
























-1 


30 


7922.464 


1.54E-04 



























36 


7928.100 


3.25E-04 
























-1 


32 


7946.406 


1.54E-04 
























-1 


34 


7971.783 


1.54E-04 
























-1 


36 


7998.571 


1.54E-04 
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Table A.3.a CN A^n - X^S+ Red System (2,0) band. 

B T' W A (air) [A] fjJ'J") EW [mA] 

*080 *102 *235 



6 


8.5 


9 


7952.261 


4.49E-05 


6 


7.5 


8 


7946.181 


4.65E-05 


6 


6.5 


7 


7940.433 


4.71E-05 


6 


5.5 


6 


7935.020 


4.69E-05 


3 


9.5 


9 


7934.855 


1.36E-04 


5 


8.5 


9 


7934.812 


1.22E-04 


3 


8.5 


8 


7930.849 


1.27E-04 


5 


7.5 


8 


7930.811 


1.30E-04 


6 


4.5 


5 


7929.943 


4.46E-05 


3 


7.5 


7 


7927.167 


1.18E-04 


5 


6.5 


7 


7927.134 


1.40E-04 


6 


3.5 


4 


7925.204 


3.84E-05 


3 


6.5 


6 


7923.813 


1.08E-04 


5 


5.5 


6 


7923.784 


1.51E-04 


6 


2.5 


3 


7920.804 


2.56E-05 


3 


5.5 


5 


7920.786 


9.53E-05 


5 


4.5 


5 


7920.762 


1.61E-04 


3 


4.5 


4 


7918.091 


8.04E-05 


5 


3.5 


4 


7918.071 


1.70E-04 


3 


3.5 


3 


7915.729 


6.22E-05 


5 


2.5 


3 


7915.714 


1.75E-04 


2 


9.5 


9 


7915.407 


3.49E-04 


4 


8.5 


9 


7915.364 


7.82E-05 


3 


2.5 


2 


7913.704 


3.64E-05 


5 


1.5 


2 


7913.692 


1.64E-04 


2 


8.5 


8 


7913.473 


3.40E-04 


4 


7.5 


8 


7913.435 


8.68E-05 


9 


8.5 


9 


7913.336 


1.49E-04 


2 


7.5 


7 


7911.856 


3.29E-04 


4 


6.5 


7 


7911.822 


9.73E-05 


2 


6.5 


6 


7910.559 


3.18E-04 


4 


5.5 


6 


7910.529 


l.lOE-04 


2 


5.5 


5 


7909.584 


3.05E-04 


4 


4.5 


5 


7909.559 


1.26E-04 


2 


1.5 


1 


7908.966 


2.08E-04 


4 


0.5 


1 


7908.959 


4.37E-04 


2 


4.5 


4 


7908.934 


2.91E-04 


4 


3.5 


4 


7908.914 


1.48E-04 


2 


2.5 


2 


7908.622 


2.50E-04 


2 


3.5 


3 


7908.613 


2.73E-04 


4 


1.5 


2 


7908.611 


2.45E-04 


4 


2.5 


3 


7908.597 


1.81E-04 


9 


7.5 


8 


7908.172 


1.43E-04 


1 


0.5 





7906.598 


4.93E-04 


1 


1.5 


1 


7903.892 


3.12E-04 


9 


6.5 


7 


7903.237 


1.37E-04 


1 


2.5 


2 


7901.520 


2.59E-04 


1 


3.5 


3 


7899.481 


2.37E-04 


9 


5.5 


6 


7898.529 


1.29E-04 


1 


4.5 


4 


7897.771 


2.26E-04 


1 


5.5 


5 


7896.386 


2.21E-04 


1 


6.5 


6 


7895.326 


2.18E-04 



17.3 
43.9 



65.1 39.4 



109.1 52.2 



36.1 37.9 36.9 

47.2 71.5 18.2 

30.9 23.1 
33.5 

38.5 
15.7 
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Table A.S.b Continued; CN A^n - X^E+ Red System (2,0) band. 
"B T' W A (air) [A] fjJ'J") EW [mA] 

*080 *102 *235 



12 


9.5 


9 


7895.174 


6.22E-05 


8 


8.5 


9 


7895.131 


3.42E-04 


1 


7.5 


7 


7894.583 


2.17E-04 


1 


8.5 


8 


7894.158 


2.17E-04 


1 


9.5 


9 


7894.045 


2.17E-04 


9 


4.5 


5 


7894.043 


1.21E-04 


12 


8.5 


8 


7892.153 


6.68E-05 


8 


7.5 


8 


7892.116 


3.32E-04 


9 


3.5 


4 


7889.776 


l.llE-04 


12 


7.5 


7 


7889.360 


7.14E-05 


8 


6.5 


7 


7889.327 


3.22E-04 


12 


6.5 


6 


7886.792 


7.60E-05 


8 


5.5 


6 


7886.763 


3.11E-04 


9 


2.5 


3 


7885.725 


9.84E-05 


12 


5.5 


5 


7884.445 


8.06E-05 


8 


4.5 


5 


7884.420 


2.99E-04 


12 


4.5 


4 


7882.315 


8.46E-05 


8 


3.5 


4 


7882.295 


2.88E-04 


9 


1.5 


2 


7881.889 


7.79E-05 


12 


3.5 


3 


7880.400 


8.78E-05 


8 


2.5 


3 


7880.384 


2.77E-04 


12 


2.5 


2 


7878.697 


8.76E-05 


8 


1.5 


2 


7878.686 


2.70E-04 


12 


1.5 


1 


7877.205 


7.79E-05 


8 


0.5 


1 


7877.198 


3.09E-04 


11 


9.5 


9 


7874.878 


1.14E-04 


11 


0.5 





7874.852 


3.09E-04 


7 


8.5 


9 


7874.835 


1.93E-04 


11 


1.5 


1 


7873.992 


2.26E-04 


11 


8.5 


8 


7873.987 


1.23E-04 


7 


0.5 


1 


7873.985 


1.84E-04 


7 


7.5 


8 


7873.949 


1.90E-04 


11 


2.5 


2 


7873.343 


2.01E-04 


7 


1.5 


2 


7873.332 


1.73E-04 


11 


7.5 


7 


7873.325 


1.34E-04 


7 


6.5 


7 


7873.292 


1.86E-04 


11 


3.5 


3 


7872.905 


1.84E-04 


7 


2.5 


3 


7872.889 


1.73E-04 


11 


6.5 


6 


7872.889 


1.45E-04 


7 


5.5 


6 


7872.860 


1.82E-04 


11 


4.5 


4 


7872.682 


1.70E-04 


11 


5.5 


5 


7872.676 


1.57E-04 


7 


3.5 


4 


7872.662 


1.75E-04 


7 


4.5 


5 


7872.652 


1.78E-04 


10 


0.5 





7871.654 


1.28E-04 


10 


1.5 


1 


7868.667 


1.07E-04 


10 


2.5 


2 


7865.892 


9.53E-05 


10 


3.5 


3 


7863.333 


8.60E-05 



23.5 18.7 
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Table A.4.a CN A^n - X'^E+ Red System (3,0) band. 





J 


iV 


A (airj [AJ 


/(j J ) 








ii W 


A 1 

|mA] 


















*042 


*051 


*053 


*561 


*200 


*268 


*222 


*233 


6 


8.5 


9 


6987.583 


1.98E-05 


















6 


7.5 


8 


6982.767 


2.04E-05 


















6 


6.5 


7 


6978.223 


2.08E-05 


















3 


9.5 


9 


6974.278 


5.93E-05 


















5 


8.5 


9 


6974.244 


5.33E-05 


















6 


5.5 


6 


6973.953 


2.05E-05 


















3 


8.5 


8 


6971.045 


5.57E-05 


















5 


7.5 


8 


6971.016 


5.74E-05 


















6 


4.5 


5 


6969.958 


1.96E-05 


















3 


7.5 


7 


6968.079 


5.18E-05 


















5 


6.5 


7 


6968.053 


6.18E-05 


















6 


3.5 


4 


6966.241 


1.68E-05 
















6.5 


3 


6.5 


6 


6965.383 


4.71E-05 


















5 


5.5 


6 


6965.360 


6.62E-05 


















3 


5.5 


5 


6962.956 


4.18E-05 


















5 


4.5 


5 


6962.937 


7.08E-05 










79.7 








6 


2.5 


3 


6962.801 


1.12E-05 


11.5 
















3 


4.5 


4 


6960.801 


3.53E-05 


















5 


3.5 


4 


6960.785 


7.47E-05 


8.2 








70.3 




73.2 


53.2 


2 


9.5 


9 


6959.402 


-1 P^ OTTl f\ A 

1.53E-04 


















4 


8.5 


9 


6959.369 


3.45E-05 


















3 


3.5 


3 


6958.920 


2.70E-05 


















5 


2.5 


3 


6958.908 


7.68E-05 


11.5 


25.6 


46.8 


12.8 


76.6 


19.2 


72.9 


43.2 


2 


8.5 


8 


6957.753 


1.48E-04 


















4 


7.5 


8 


6957.724 


3.83E-05 


















9 


8.5 


9 


6957.622 


6.54E-05 


















3 


2.5 


2 


6957.314 


1.60E-05 


















5 


1.5 


2 


6957.305 


7.19E-05 


19.3 


42.2 




7.9 


58.6 


30.7 


67.4 




2 


7.5 


7 


6956.367 


-1 A A 1 ^ J 

1.44E-04 


















4 


6.5 


7 


6956.340 


4.28E-05 


















2 


6.5 


6 


6955.243 


1.39E-04 


















4 


5.5 


6 


6955.220 


4.84E-05 


















2 


5.5 


5 


6954.385 


1.33E-04 


















4 


4.5 


5 


6954.366 


5.55E-05 


















2 


4.5 


4 


6953.795 


1.27E-04 


















4 


3.5 


4 


6953.779 


6.50E-05 


















2 


1.5 


1 


6953.652 


9.13E-05 


















4 
9 


0.5 
7.5 


1 

8 


6953.646 
6953.498 


1.92E-04 
6.27E-05 




58.2 














2 


3.5 


3 


6953.475 


1.20E-04 


















4 


2.5 


3 


6953.462 


7.99E-05 


















2 


2.5 


2 


6953.427 


l.lOE-04 


















4 


1.5 


2 


6953.418 


1.08E-04 


















1 


0.5 





6951.821 


2.15E-04 


70.9 


106.8 




34.1 


122.9 


104.2 


145.2 


110.8 


1 


1.5 


1 


6949.770 


1.36E-04 


15.2 




34.1 


24.0 


101.0 


114.4 


128.9 


55.7 


9 


6.5 


7 


6949.568 


5.98E-05 


















1 


2.5 


2 


6947.992 


1.14E-04 


20.9 




42.0 


14.3 


152.5 


66.9 


115.7 


44.3 


1 
1 


0.0 


q 





1 r\AT? r\A 














01.1 


qv 7 


9 


5.5 


6 


6945.831 


5.66E-05 


















1 


4.5 


4 


6945.251 


9.91E-05 


















1 


5.5 


5 


6944.285 


9.66E-05 






67.0 












12 


9.5 


9 


6943.727 


2.74E-05 
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Table A.4.b Continued; CN A^n - X^T,+ Red System (3,0) band. 





J" 
J 


\T" 
iV 


A (airj [AJ 


ikJ J ) 








771T17 r. X 1 

rjW [mAJ 


















*042 


*051 


*053 


*ODl *zOO 


*268 


*222 


*233 


8 


8.5 


9 


6943.693 


1.50E-04 
















1 


6.5 


6 


6943.584 


9.53E-05 
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Table A.5. CH+ A^U - X^T,+ (0,0) band. 
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